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We consider charge qubits based on shallow donor electron states in silicon and coupled quantum 
dots in GaAs. Specifically, we study the feasibility of P J charge qubits in Si, focusing on single qubit 
properties in terms of tunnel coupling between the two phosphorus donors and qubit decoherence 
■ caused by electron-phonon interaction. By taking into consideration the multi-valley structure of 

(~| ' the Si conduction band, we show that inter-valley quantum interference has important consequences 

for single-qubit operations of charge qubits. In particular, the valley interference leads to a 
tunnel-coupling strength distribution centered around zero. On the other hand, we find that the 
Si bandstructure does not dramatically affect the electron-phonon coupling and consequently, qubit 
coherence. We also critically compare charge qubit properties for Si:Pj and GaAs double quantum 
dot quantum computer architectures. 



C/3 



o 



> 



o 



■ 

X 



PACS numbers: 71.55.Cn, 03.67.Lx, 85.35.-p 



I. INTRODUCTION 



Among the solid state candidates for qubits in quantum information processing, semiconductor-based systems have 
been among the most extensively explored. Key features in favor of these proposals are the high level of theoreti- 
cal understanding, experimental control, and nanofabrication capabilities currently available for semiconductors. It 
i is commonly believed that group-IV or III-V semiconductor nanostructure-based quantum computer architectures 
' should be relatively easily scalable because of the existence of the vast semiconductor microelectronics infrastructure. 
: This scalability incentive has led to a great deal of recent activities in studying qubit properties of semiconductor 
nanostructures.^'^ Theoretically, many semiconductor-based quantum computers were conceived to rely on either 
electron spins or nuclear spins as qubits li^i^iSii^iSiifl Spin-1/2 fermions (electrons or nuclei) probably constitute the 
most natural and robust choices of quantum two-level systems for qubits in solids. Unfortunately, and in spite of 
considerable recent progress, "'^■'^'■'^^ it is difficult to perform fast measurements of single electron (or nuclear) spins, 
. which are required for practical quantum computing implementations incoporating quantum error corrections. The 
problem here is quantitative-although the electron spin is surely a quantum two-level system, the Bohr magneton is a 
i' , very small quantity, (and spin usually does not couple strongly to external probes), making it difhcult to manipulate 
'■^ and measure the microscopic single spin states in the solid state environment. 
^ [ In contrast to spin qubits, charge qubits in semiconductors have the substantial advantage of being easy to ma- 
O • nipulate and measure since the experimental techniques for measuring single electron charges in semiconductors are 
] extremely well-developed. The price one pays for the relative ease in the manipulation and read-out of single-charge 
K*" • states is, of course, the strong decoherence and the rather short decoherence time of the orbital charge states be- 
cause they couple strongly to the environment through the long-range Coulomb interaction. This fast decoherence 
of orbital states makes semiconductor charge qubits rather unlikely candidates for a scalable quantum computer ar- 
il^ ■ chitecture. However, the strong interactions make the orbital states an excellent choice for studying qubit dynamics 
and qubit coupling in the solid state nanostructure environment. This is particularly true in view of the difficulties 
encountered in the manipulation and the measurement of the single spin states in semiconductors. It is worthwhile 
also to remember that the much-studied superconducting-Cooper-pair-box-based quantum computer architectures 
are charge-based systems as well^^ii^ and there are conceptual and formal overlaps between semiconductor charge 
qubits and superconductor charge qubits, providing further impetus for studying orbital qubits in semiconductor 
nanostructures. 

There have been several proposals for orbital/charge qubits in semiconductorsii^4iS*iii2ii2i2£ In this work we the- 
oretically analyze single charge qubit properties for phosphorus donor states in silicon, comparing it critically with 
charge qubit states in coupled quantum dots (QD) in GaAs. Our specific goal is to investigate how the peculiar six-fold 
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valley degeneracy of Si conduction band affects the single qubit properties of orbital states in P-doped Si system. 
The issue is important in the context of our earlier results showing that quantum interference between valleys leads 
to a strong suppression of the exchange energy controlling the inter-qubit coupling in the electron-spin-based silicon 
quantum computer architecturei^i*^ The dramatic adverse effect of the valley interference on the silicon exchange 
gate naturally raises the question of whether a similar valley interference effect would also strongly (and adversely) 
affect the charge qubit properties. We answer this question in this paper. 

Successful coherent manipulation of electron orbital states in GaAs has been achieved for electrons bound to donor 
impuritiesS2i as well as electrons in double quantum dots.— ^ There were also suggestions of directly using electron 
orbital states in GaAs or Si as the building blocks for quantum information processingiiSiiLi2i2& Specifically, a double 
QD with an electron bound in each dot or a pair of phosphorus donors that sit relatively close to each other (so as to 
have sizable wave function overlap) form an effective hydrogen molecule in GaAs or Si host material. Charge qubits 
may be defined by ionizing one of the bound electrons, thus leading to a double well potential filled with a single 
electron: The single electron ground state manifold, whether it is the two states localized in each of the wells or their 
symmetric and anti-symmetric combinations, can then be used as the two-level system forming a charge qubit»iS*ii 
The advantage of such a charge qubit is that it is easy to manipulate and detect, while its disadvantage, as already 
mentioned above, is the generally fast charge decoherence as compared to spin. 

Here we study the feasibility of the charge qubit in Si, focusing on single qubit properties in terms of the tunnel 
coupling between the two phosphorus donors (Sec^, and charge decoherence of this system in terms of electron- 
phonon coupling fSec lIII|l . We take into consideration the multi- valley structure of the Si conduction band and explore 
whether valley interference could lead to potential problems or advantages with the operations of Pj charge qubits, 
such as difhculties in the control of tunnel coupling similar to the control of exchange in two-electron systemsj^LSS 
or favorable decoherence properties through vanishing electron-phonon coupling. In section Hvl we critically compare 
charge qubits based on Si:Pj and GaAs double QD systems. 



II. THE SYMMETRIC-ANTISYMMETRIC GAP FOR THE P+ MOLECULE IN SILICON: QUBIT 

FIDELITY 

We study the simple situation where a single electron is shared by a donor pair, constituting a Pj molecule in 
Si. The charge qubit here consists of the two lowest energy orbital states of an ionized P2 molecule in Si with 
only one valence electron in the outermost shell shared by the two P atoms. The key issue to be examined is the 
tunnel coupling and the resulting coherent superposition of one-electron states, rather than the entanglement among 
electrons, as occurs for an exchange-coupled pair of electrons. 

The donors are at substitutional sites R^i and Rb in an otherwise perfect Si structure. In the absence of an 
external bias, and for well separated donors, we may write the eigenstates for the two lowest-energy states as a 
superposition of single-donor ground state wavefunctions localized at each donor, ipAi'r) and V'-B(r), similar to the 
standard approximation for the molecular ioniS The symmetry of the molecule leads to two eigenstates on this 
basis, namely the symmetric and antisymmetric superpositions 

v/2(l±5) ^ ' 

where S = is the overlap integral and is real.^^ The conduction band of bulk Si has six degenerate minima 

{n = 1,...,6), located along the F— X axis of the Brillouin zone at |k^| ~ 0.85(27r/a) from the F point, where 
a = 5.43 A is the Si lattice parameter. Following Kohn-Luttinger effective mass approximation)S& the single-donor 
ground state wavefunctions are written in terms of the six unperturbed Si band edge Bloch states (f>^ = u^(r)e''''' ''. 
For the donor at Ra, 



^A{r) = Y,^M^-^A)Mr,-RA) = ^«^i^^(r-R^)u^(r)e'''"-('-^^) , (2) 

where the envelope functions centered at the donor site, -F)i(r — Ra), are deformed shallow donor effective mass 
IS hydrogenic orbitals. For instance, for ^ = z, -F'z(r) = exp{ — [(x^ -I- y^^jc? -I- j}?^!'^^ j \J i^a^h. The expansion 
coefficients a^, which are also called valley populations, are realiSi The effective Bohr radii a and h are variational 
parameters chosen to minimize Ea = (^Ra I-^^^IV'Ra): leading to a = 25 A and 6 = 14 A when recently measured 
effective mass values are used in the minimization^ The operator Ha (note that in our notation the single donor 
Hamiltonians Ha and Hb are equivalent) is the single-donor Hamiltonian for the bound electron^ which includes the 
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FIG. 1: (color online) Symmetric-antisymmetric gap for the molecular ion in Si for the donor pair along the indicated 
lattice directions. The arrow in the upper frame indicates the target configuration analyzed in Fig. |5| 



kinetic energy, the Si periodic potential, the impurity screened Coulomb potential centered at Ra, and the valley-orbit 
effects, leading to Ea ~ —40 mcV, consistent with the experimentally observed value of 45 meV. 
The Hamiltonian for the singly ionized donor pair can then be written as 

2 2 

from which it is straightforward to obtain the expectation value of the single qubit energy gap between the lowest-lying 
(symmetric and antisymmetric) states 4'±(r) in 

2 ® 

As-AS=<*-|i^|*- >"<*+|i?|*+>= J— ^^<5^(R)cos(k^-R) , (4) 

where = l/\/6 for unstrained Si;^ i^/i(R-) \c(u\'^(suCi — Vfj.) and expressions for S, Ci, s^j and Vf^, all of which are 
functions of R = Ra — Rb, are given in Ref. |22. For |R| ^ a, 6, 5* <C 1, and the amplitudes ^^(R) are monotonically 
decaying functions of the interdonor distance R. Except for the anisotropy, which is a consequence of the effective 
mass anisotropy in Si, the dependence of 5^ on |R| is qualitatively similar to the symmetric-antisymmetric gap in 
the molecule, namely an exponential decay with power-law prefactors. The main difference here comes from the 
cosine factors, which are related to the oscillatory behavior-^ of the donor wavefunction in Si arising from the Si 
conduction band valley degeneracy and to the presence of two pinning centers. 

Figure^gives the calculated gaps as a function of R for a donor pair along three high-symmetry crystal directions. 
Two points are worth emphasizing here, which are manifestly different from the corresponding hydrogenic molecular 
ion behavior: (i) As_as is an anisotropic and fast oscillatory function of R; (ii) the sign of A5_^5 may be positive or 
negative depending on the precise value of R. The characteristics mentioned in point (i) are similar to the exchange 
coupling behavior previously discussed for the two-electrons neutral donor pairiSLSStS^ Point (ii) implies that the 
molecular ion ground state in Si may be symmetric (as in the molecular ion case) or antisymmetric depending on 
the separation between the two P atoms. Note that for the two-electron case, the ground state is always a singlet (i.e. 
a symmetric two-particle spatial part of the wavefunction with the spin part being antisymmetric), implying that the 
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FIG. 2: (color online) Probability distributions of the symmetric-antisymmetric gap for the molecular ion in Si. Donor 
pairs are approximately aligned along [100] , but with an uncertainty radius Ru with respect to this target axial alignement (see 
text). The arrow in each panel indicates the gap value for the target configuration, for which the uncertainty radius is Ru ~ 0. 
Notice that all four distributions are peaked at As_as = 0, and not at the target gap value. 



exchange J is always positive for a two-electron molecule. For a one-electron ionized molecule, however, the ground 
state spatial wavefunction can be either symmetric or antisymmetric. 

Figure 13 shows normalized probability distributions for the Ag-^s gap values when the first donor is kept fixed 
at Ra and the second donor is placed at a site 20 lattice parameters away (~ 108.6 A), along the [100] axis. This 
target configuration is indicated by an arrow in Fig. We allow the second donor position Rs to visit all possible 
substitutional diamond lattice positions within a sphere of radius Ru centered at the attempted position. Our 
motivation here is to simulate the realistic fabrication of a molecular ion with fixed inter-atomic distance in Si 
with the state of the art Si technology, in which there will always be a small {Ru ~ 1 — 3 nm) uncertainty in the precise 
positioning of the substitutional donor atom within the Si unit cell. We would like to estimate the resultant randomness 
or uncertainty in Ag-^s arising from this uncertainty in R^. For i?„ — 0, i.e., for R = 20ai, Ag^^is ~ 2.4 meV, 
given by the arrows in Fig. |21 We incorporate the effect of small uncertainties by taking i?„ = Inm, corresponding to 
the best reported degree of accuracy in single P atom positioning in Si3S These small deviations completely change 
the qubit gap distribution, as given by the histogram in Fig.|2Ia) strongly peaked around zero. A similar distribution 
is obtained for i?„ — 2 nm, as illustrated in Fig.jS^b). Further increasing Ru leads to broader distributions of the gap 
values, though still peaked at zero [see Fig.j^^c) and (d)]. This broadening is due to the fast increase in the number of 
lattice sites inside the sphere of radius Ru, thus contributing to the distribution, as Ru increases. We conclude that 
the valley interference between the six Bloch states leads to a strong suppression of the qubit fidelity since the most 
probable Ag-^g tends to be zero. 

A very small As-as is undesirable in defining the two states |0) and |1) forming the charge qubit. If we take 
them to be the symmetric and anti-symmetric states given in Eq. the fact that they are essentially degenerate 
means that, when one attempts to initialize the qubit state at |0), a different combination a|0) -|- might result. 
Well defined qubits may still be defined under a suitable applied external bias, so that the electron ground state 
wavefunction is localized around one of the donors, say at lattice site R^, and the first excited state is localized arond 
Rs- 

Single qubit rotations, used to implement universal quantum gates, might in principle be achieved by adiabatic 
tunneling of the electron among the two sites under controlled axially aligned electric fields through bias sweepsi^ 
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When, at zero bias, the ground state is not well separated by a gap from the first excited state, severe limitations are 
expected in the adiabatic manipulation of the electron by applied external fields. In other words, the fidelity of the 
single qubit system defining the quantum two-level dynamics will be severely compromised by the valley interference 
effect. 



III. ELECTRON-PHONON COUPLING FOR A P+ MOLECULE 

Two key decoherence channels for charge qubits in solids are background charge fluctuations and electron-phonon 
coupling.^'* The former is closely related to the sample quality (e.g., existences of stray charges and charged defects 
in the system) and is extrinsic, while the latter is intrinsic. Here we focus on the electron-phonon coupling. A critical 
question for the P J molecular ion in Si is whether the Si bandstructure and the associated charge density oscillations^^ 
lead to any significant modification of the electron-phonon coupling matrix elements. In the following, we first derive 
the electron-phonon coupling for a single valley situation, such as GaAs, to set a benchmark, then assess the effect of 
the Si conduction band valleys and the Bloch functions on the donor electron-phonon coupling matrix elements. Our 
motivation is to investigate whether valley interference leads to a strong suppression of the electron-phonon coupling 
(similar to the suppression of exchange and tunnel couplings), which would be beneficial for silicon charge qubits. 

According to the results in Sec^ the energy splitting Ag-^s between the two lowest energy states in a molecular 
ion is up to a few meV, thus only low energy acoustic phonons near the Brillouin zone center q ~ contribute to 
phonon-induced electron decoherence. Electron-acoustic-phonon coupling in a semiconductor can be classified into 
two types: deformation potential and piezoelectric. Since Si is not polar, deformation potential is the only relevant 
interaction. We will thus focus on this interaction for the rest of this paper. The eletron-phonon interaction takes the 
form:^ 

He,-DY^ {^^) |q|p(q)(«. + «!,) , (5) 

where D is the deformation constant, p,„ is the mass density of the host material, V is the volume of the sample, Oq 
and are phonon annihilation and creation operators, and p(q) is the Fourier transform of the electron density 
operator: 

p(q) = ^ 4 c, I dr e-^^-r, (r)0, (r) , (6) 

A, 17 

where A and 77 are indices of electronic states/modes, c\ and c\ are electronic annihilation and creation operators for 
the A-mode, while (j) are mode functions. For our two-donor (or double-dot) situation, where we are only interested in 
the two lowest energy single-electron eigenstates, we can choose them as the basis (so that A and 77 take the value of 
+ and — as defined below) and the electron-phonon coupling Hamiltonian is conveniently written in this quasi-two- 
level basis in terms of the Pauli spin matrices ax and az (where spin up and down states refer to the two electronic 
eigenstates): 

/ h \^/^ 



q 



A^ - ^((+|e^''--|+>-(-|e^''--|-» . (7) 

Here the term proportional to ax can lead to transition between the two electronic eigenstates and is related to 
relaxation; while the term proportional to ct^ only causes energy renormalization of the two electronic levels, but no 
state mixing, so that it only leads to pure dephasing for the electronic charge states. 

We first consider a double quantum dot with one electron (which is similar to a singly ionized donor pair) in GaAs 
where the nondegenerate conduction band minimum occurs at the F point. When the two dots or donors are well 
separated and not strongly biased, the relevant single-electron states are 

$+ = a(j)A{r) + b(j>B{r) ; = ^^^(r) - a0B(r), (8) 

with 0A(B)(r) — </'(r — Ryi(B)) ■"o(i")j where (p{'r) is a slowly varying envelope function, and the Bloch function at 
the conduction band minimum (k = at F point) is equal to the periodic part uo{r). Though we have chosen the 
envelopes if centered at each well to be identical, they could be different, as is generally the case for quantum dots. 
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The deformation potential electron-phonon coupling matrix element for relaxation between the initial unperturbed 
eigenstate |— ) and the final eigenstate |+) is proportional to (— |e*'i '"|+), where q is the phonon wavevector. For a 
GaAs double donor or double dot case the matrix element is proportional to 

Ar - ($_|e'i"-|$+) = j l^oWpe*'!-'' {ab*[ip{r)f - a*b[ip{Y - B)f 

+ (|fep-|apMrMr-R)} . (9) 

For small energy splittings between the $± states, all terms in the integrand of Eq. Q are slowly varying functions in 
the interatomic spacing scale, except |uo(r)p, which is periodic and normalized in a primitive cell of the host material: 
^ (ir|uo(r)P = 1 where Vt is the volume of the primitive cell. This allows for the approximation J dr |ti(r)p/(r) « 
/ dr f{r),^^ valid for slowly varying /(r), to be applied to ||2J), leading to: 

Ar = {ab* -a*be''^ '^) J dr e''i '"[</j(r)]2 

+ - lap) J dr e''i-Xr)(p(r - R) . (10) 

Here the first integral is an on-site contribution modified by the phase difference e"i'^ between the two dots/donors, 
while the second integral is a two-dot contribution that is generally much smaller because of the small overlap. 
The dephasing matrix element A^ can be similarly calculated and the result is 

A^ = z(|&|2_|a|2)eW2sin3_5: j dr e^^ ^^^ir)]^ 

+ {a*b + ab*) j dr e'"^''' ip{r)ip{r - R) . (11) 

Notice that here the prefactors |&p — |ap and a*b + ab* are for on-site and off-site integrals, just the opposite to 
what we have in Eq. (|10|) . In other words, when |6| ~ |a|, A^f, is small, charge decoherence caused by electron-phonon 
interaction is dominated by relaxation;'^^'"^'* when \b\ and \a\ are very different (so that, for example, |6| ~ 1 and 
\a\ ^ 0), charge decoherence is dominated by pure dephasing^ Below we will focus on the relaxation matrix element 
Ar = (— |e*'' '"|-|-) as the contributing integrals are identical in the dephasing matrix element A^. 

We now consider a singly ionized phosphorus donor pair in Si, taking into account the Si bandstructure. For 
two donors not too close to each other, and possibly detuned by an axially aligned electric field, the lowest energy 
single-electron states are superpositions of V'A centered at Ra [given in Eq. (j^J] and ijjB centered at Rs, similar to 
Eq. ®: 

*+ = aijAir) + bi^Bir) ; = bi^A{r) - aipsir), (12) 

where the superposition coefficients a and b are of course not to be confused with the effective Bohr radii. If we choose 
R^ = 0, Rb — R, and as the initial and final states, the relaxation matrix element Ar can be written as 

(M'_|e^i-'-|*+) = dr u*^{r)u,{r)e-'^''--''-^-^>'' {ab*Ff,{r)F,{r) 

-a*bF^{r - R)i^^(r - R)e»(kM-k.)-R + \b\^F^{r)F^{r - R)e-*''-^ 
-|a|2F,(r)F^(r~R)e'''--^} . (13) 

If the energy splitting for the two double donor states is small (because of large inter-donor separations or valley 
interference), the dominant electron-phonon coupling is restricted to the regime of |q| ^ |k^|. A simpler form 
of Eq. (|13|l can then be obtained by just keeping those integrals with ji = v (other integrals have fast oscillatory 
integrands and are thus vanishingly small): 

6 

(*_|e*i '"|*+) = {ab* - a*b e*i-^) ^ ja^p / dr e^'^ '' Fl{r) 

6 

H\b? ~ l«n E l«A'l'cos(k^ • R) y dr e'^-''F^{r)F^{r - R) . (14) 
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Equation 1)14(1 for double donor state in Si takes on quite similar form as Eq. H10|l for double dot states in GaAs, 
with the first sum in Eq. 114|) containing on-site contributions, and the second sum containing off-site (inter-dot) 
contributions. Thus the first sum should generally outweigh the second even for a ^ b. However, if one of the 
coefficients a or 6 is very small, for example due to electric field bias, the second sum may become dominant for the 
relaxation matrix element. However, as we mentioned above, in that situation pure dephasing becomes the more 
important decoherence channel. In the case when the overlap integrals do make non-negligible contributions (for 
example, when the two donors are detuned but not too strongly so, and the two donors are sufficiently close so that 
the overlap integrals are not vanishingly small), it is interesting to note that each of the terms in the sum over the 
valleys is multiplied by the same cos(k^ • R) factors which appear in Eq. Q. The effect again is to have results for 
the off-site contribution to the electron-phonon coupling strongly oscillatory as a function of the interdonor relative 
position R. The average overall effect, as illustrated in Fig. |21 is to reduce the absolute value of the relaxation 
coupling. 

We now consider the possible contributions when |q| may not be negligibly small. Indeed, in Si, for hcuq ~ 5 meV, 
q - 0.1^. Thus, if the energy splitting between 'if± states is ^ 5 meV, we need to include in our calculation phonon 
wave vectors that may couple to the periodic part of the Bloch functions as described below.^^ Expanding the periodic 
part of the Bloch functions it^ in terms of plane waves (restricted to the reciprocal lattice wave vectors) yields: 



iG„r 



SO that the relaxation matrix elements of Eq. (|13|l becomes 

X [ab*F^{v)F,{v) ~ a*hF^{v - R)F,(r - R)e^(k.-k.)-R 

~\a\^F^{v ~ K)F,{vy--^ + \h\^F^{v)F,{v - R)e-^''-^] . (15) 

Since |q| is always relatively close to zone center (i.e., |q| is always much smaller than 27r/a), the largest contribution 
to Eq. ((15|l comes from terms with v — ^, and = . These are the same terms that determine the matrix elements 
in the small q limit, as given by Eq. (|14|) . The important question now is whether other terms will also contribute 
significantly when |q| is not particularly close to the zone center. Recall that more than 90% of the spectral weight 
in comes from five plane waves^* (the rest of the Cq coefficients are at least one order of magnitude smaller): 
For u^, these are G^; = 0, ^(-1,±1,±1), so that K + G^k. 0.85 ^, ^(-0.15, ±1, ±1) are the five smallest wave 
vectors contributing to the Bloch function ii2,(r)e''^^''=^ . There is thus one scenario when -I- k^ — K^, — ki, — q might 
have similar amplitude as q: when v = and G^ — K,y is parallel to k^. For example, for q = qx, there are terms 
with Gx + kj, — K_2, — k_j. — ^(±0.3, 0, 0). Since these wave vectors correspond to wave lengths of the order of 15 
A, while the donor effective Bohr radius is about 20 A, one needs to carefully evaluate the integrals involving these 
terms as their oscillatory integrands have the same length scale as the envelopes. The = —fi contribution to the 
electron-phonon coupling matrix elements takes the form (taking into consideration that = F^, k_^ = — k^, and 

6 

(vI/_|e^'i-"-|*+).=-, = EI"mP E ^gC^-. 

X |[a6* - a*b e-*(GM-K-,,-q)-R] 

X J dr e-'(G,.+k,,-K-,,-k-,,-q).r^2(j.) _^ (|^|2 _ \a\^yK-R 

X y dr e-^iG^.+K-K-^-k-^-^yrp^i^j..^p^(^j. _ R)| . (16) 

For each /i, only the 5 dominant G^ contributions mentioned previously are included in the second summation above. 
(As an example, in Table |l] we list dataS^ for /i — ±x). For each G^ ^ 0, there exists one and only one K_^ for 
which the exponent G^^ -I- k^ — — k_^ is small. If G^ = 0, all exponents are large so that the integrals should 
not be important in the sum in Eq. H16|l. Without loss of generality, we consider the fji = x part of the sum, where 
Gx + kj; — K_j; — k_2, — ^(—0.3,0,0). Notice that this exponent is independent of which G^; is in consideration. 
Therefore, all the terms except Cq C'k_j, can be factored out of the sum over G^ and K_^, so that the sum in 
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Eq. ((TB|) is proportional to X^g^ k^^^o ^g^^^-^^ which vanishes due to the symmetry of Si lattice, as illustrated in 
TableQ] Therefore, for intermediate |q| the lowest order correction to the small |q| electron-phonon coupling matrix 
element vanishes, so that Eq. 114(1 is valid in both small and intermediate |q| regimes for the phonons involved. 





(-1,1,1) (-1,1,-1) (-1,-1,1) (-1,-1,-1) 
(0.343,0) (-0.313, +0.313) (-0.313, -0.313) (-0.313, -0.313) (-0.313, 0.313) 




(1,1,1) (1,1,-1) (1,-1,1) (1,-1,-1) 
(0.343,0) (-0.313, -0.313) (-0.313, 0.313) (-0.313, 0.313) (-0.313, -0.313) 




(0.118,0) (0,0.196) (0,-0.196) (0,-0.196) (0,0.196) 



TABLE I: Here we give the 5 most important expansion coefficients for the periodic part of the Bloch state (px and (/>_i [for 
example, 4ix{r) = (7j-,^g»(Gx-i-ka,) rj^ Notice that C'a^ and Ck_^ are complex in general. 

In summary, the electron-phonon coupling for a molecular ion in Si formally behaves very similarly to that 
for a single electron trapped in a GaAs double quantum dot (restricted to the deformation interaction). The more 
complicated multi- valley bandstructure of Si and the strong inter- valley coupling introduced by the phosphorus donor 
atoms do not cause significant changes in the electron-phonon coupling matrix elements. The only valley interference 
effect occurs when the overlap between the two donors is not negligible. Even then the interference among the valleys 
only causes oscillatory suppression of the off-site contributions, which are relatively small anyway. Therefore, available 
estimates^i*^ of decoherence induced by electron-phonon coupling based on a single- valley hydrogenic approximation 
in the system in Si should be valid. In other words, the multi-valley quantum interference effect does not provide 
any particular advantage (or disadvantage) for single qubit decoherence in the Si:P donor charge-based quantum 
computer architecture. 

IV. DISCUSSION AND SUMMARY 

We have so far explored the feasibility of charge qubits based on the Pj system in Si. We find that this system 
possesses decoherence properties (induced by electron-phonon coupling) similar to a GaAs double quantum dot. The Si 
bandstructure does, however, significantly (and adversely) influence the tunnel coupling between the two phosphorus 
donors, so that for many relative positions between the donors in a pair, the tunnel coupling becomes quite small. 
In other words, if two donors are randomly placed in a Si host, keeping their distance approximately constant, the 
tunnel coupling between the donor sites can vary over a wide range of values (peaked around zero) because of the Si 
conduction band valley degeneracy. This is obviously rather bad news for charge qubits in the P^ system in Si: It 
implies that a large percentage of the fabricated charge qubits are unlikely to work properly since energy splitting in 
these two-level systems is essential for quantum computation. 

The quasi-randomness of tunnel coupling in a P^ donor molecular ion in Si can be contrasted with the corresponding 
tunnel coupling in a double QD in GaAs (or Si) . The Coulomb potential of a donor provides a natural strongly localized 
confinement potential. Thus donors are really identical quantum dots, with fixed positions given by the donor nuclei, 
a fixed effective Bohr radius, and a fixed ground state energy level relative to the conduction band edge. All donor 
qubits are therefore expected to have identical properties except for the donor positioning problem. The main problem 
with donors in Si is that they break the local translational symmetry and introduce a strong valley-orbit coupling. 
Donor electron states are therefore superpositions of Bloch states from all the conduction band edges. The valley 
interference effects are thus strong in a donor system such as Si:P, leading to the atomic scale oscillations in two- 
electron exchange studied beforeSl*^ and single-electron tunneling studied here. On the other hand, a gated QD is a 
truly artificial atom, whose position, shape, size, and energy levels are all determined by the applied gate potentials 
on the metallic gates some 100 nm away. The confinement produced by gate potentials are generally quite smooth 
and shallow, and the barriers between potential minima quite broad. These slowly varying features of gated QDs 
dictate that quantum dot electronic properties are very sensitive to the tuning of applied gate potentials. It is also 
inevitable that two QDs are never identical even after careful calibration. 

As a simple illustration of the effect of uneven dots, we present in Figl^the dependence of tunnel coupling between 
a double QD as a function of inter-dot distance. We find that a 5% dot size variation leads to a 10 to 20% difference 
in the tunnel coupling. Note that, in addition to size variation, there will be inevitable fluctuations in inter-dot 
separations and barrier heights as well, leading to qubit fluctuations. Obviously, careful calibration is imperative for 
a QD system to work as a reliable charge qubit. In contrast, for donor-based charge qubits, one does not need to be 
concerned with such dot size variation problems since all P donors are identical. 

In short, both donors in Si and gated QDs in either GaAs or Si pose difficult challenges to solid state quantum 
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FIG. 3: (color online) Energy gap between the ground and first excited states of a single electron double dot as a function of 
inter-dot distance. The crosses are for two identical dots with a fixed Gaussian confinemenl!— (with a ground orbital radius 
of 12.6 nm); the squares are for situations where one dot is 5% larger; the triangles are for situations where that same dot 
is 5% smaller. At larger inter-dot distances the different-dot configurations have larger energy splittings because the dot size 
difference introduces an energy level detuning that is larger than the tunnel coupling. The inset plots the same data in the 
range of 30 to 34 nm inter-dot distance in the linear scale, where we find that the 5% dot size variation leads to a 10 to 20% 
difference in the tunnel coupling. 

information processing. For donors the challenge lies more on the fabrication process, while for gated dots the 
challenge lies more on the gating control. Which type of electron confinement (carefully calibrated gated dots or 
carefully positioned donors) may turn out to be better suited for quantum computing will ultimately be determined 
by experimental work. 

For electron decoherence we have so far limited ourselves to electron-phonon coupling. As we mentioned before, 
for charge degrees of freedom another important source of decoherence is the fluctuation in charge traps close to a 
charge qubitiiSiSi^LSS Since charge fluctuation noise can be treated in a very similar fashion as the electron-phonon 
coupling, we do not anticipate any significant qualitative difference between the Si:P system and GaAs quantum 
dots. Furthermore, since electron-phonon coupling is intrinsic, the consequent decoherence is the limit that cannot 
be improved by having better materials and fabrication quality. 

In conclusion, we show that the inter-valley quantum interference leads to a strong suppression of qubit fidelity 
in charge qubits in Si, as small nanometer-scale fluctuations in donor positioning within the Si unit cell produce 
an essentially random distribution (peaked around zero) in the energy separation between the two levels defining 
the charge qubit. We find decoherence properties of charge qubits to be qualitatively unaffected by multi- valley 
effects. For QD-based charge qubits, we find variations in qubit properties arising from fluctuations in dot sizes and 
separations, which will have to be carefully characterized individually. 
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